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Abstract 

We take the viewpoint that the physically acceptable solutions of the Lorentz- 
Dirac equation for radiation back-reaction are actually determined by a second 
i order equation of motion. We propose three methods for determining this equa- 

I tion of motion. For two example systems, the equation is obtained exactly in the 

nonrelativistic regime via each of the three approaches, and the three results are 
. found to coincide. We reveal that, for both systems, back-reaction induces a damp- 

I ing proportional to velocity and, in addition, it decreases the effect of the external 

force. 

^ ■ 1 Introduction 

The radiation back-reaction (self-force) of a point charge e having the special relativistic 
world line function r is deduced to be ([1] Ch. 16.3, [2] Ch. III.3, [3], [4]) 



where indices run from to 3, is the spacetime metric, c is the speed of light, t] = 
(2/3)e^/c^, and overdot denotes differentiation with respect to proper time (note [5]). 
This force is added to the external force /, which may depend on both the spacetime 
location and the four-velocity of the particle, / = /(r, r), to obtain the Lorentz-Dirac 
equation 

mr^ = f{r, r) + V Wk " ^) r'" (2) 



for the motion of the point particle with mass m. 



The problems with this equation are well-known. First, it is of third order so the 
initial values of spacetime position, velocity and acceleration are necessary to obtain the 
motion, and there is no physically acceptable reasoning how to prescribe acceleration. 
Second, the equation admits 'runaway' solutions — motions accelerating exponentially 
in time — , and, third, it exhibits acausal behavior. 

There are a number of attempts to treat these problems (see, e.g., [6, 7, 8, 9, 10, 11, 
12, 13, 14, 15, 16]). One of the proposed and most frequently applied way out is the 
following [17]. As a zeroth approximation, equation (2) is considered without radiation: 

mP — /-' (r, r). (3) 

The third derivative is computed from this equation: 

m \dr^ df^ J 

Then, the second derivative here is replaced by the rhs of (3) (divided by m), an expression 
of lower order derivatives, leading to the following approximation of the self-force (1): 

This is substituted into (2) to derive an approximate second order equation of motion: 

^jij — p (^j.^ f.^ _|_ jjJ^^ (^j.^ _ (•g^ 

In known examples, b^^y (r, f) is found to be a damping, dissipative term. 

Another idea is that the initial values of spacetime position, velocity and acceleration 
cannot be given independently, and one has to find a 'critical manifold' formed by those 
initial values which do not result in runaway solutions [9]. 

To see the situation more clearly, let us sharply distinguish the notion of equality 
from the notion of equation, both usually denoted by the same symbol = . An equality is 
an assertion while an equation is a definition. For instance, the equality sin^ + cos^ = 1 
asserts a relation between the functions sin and cos. The equation + 3x — l = defines 
the set of real numbers x for which this relation holds. 

An equation of motion must be an equation which defines the set of all possible 
motions under given external circumstances. 

Formula (1) is an equality, making an assertion on the self-force, for a given m,otion. 
Accordingly, we accept that the world line function r of a charged particle in an exterior 
field satisfies equality (2). This, however, does not mean that (2) is the equation of 
motion of the particle: maybe many of its solutions are not possible motions. Indeed, 
the runaway solutions are not. 

The possible motions of a radiating charged particle form a subset of the solutions of 
(2). The task is to find this subset. According to our viewpoint taken here, the set of 
possible motions must be determined by a second order differential equation. 

In what follows, we present three methods to accomplish this task. For demonstration, 
we investigate these three approaches in quantitative detail, in the nonrelativistic approx- 
imation of the Lorentz-Dirac equality, i.e., taking only j = 1, 2, 3 in (2) and omitting the 
term containing 1/c^, for two specific external force functions. 



(4) 



2 Three proposed methods 



For convenience, we introduce the shorthands 



/:=-/, V-=-V, b:^-b. (7) 
m m m 



2.1 Iteration of the radiation term 

Definitely, equation (6) can only be considered a first approximation, rather than an exact 
equation, i.e., an equation of motion whose solutions would exactly satisfy (2). In fact, if 
(6) were an exact equation then the self-force computed from its third derivative would 
give the self-force included in it. It is not hard to see that this is not so, e.g., for an elec- 
tromagnetic Lorentz force f^{r,r) = k{r)f^ in case when the external electromagnetic 
field strength F is actually constant. 

Then it is a straightforward idea that we take an analogous second approximation: 'r 
is computed as the derivative of 

= f^{r,f) + h^-y^{r,f), (8) 

and then f is substituted by the rhs of (8). Thus we obtain an expression 6(r, '^){2}' with 
which the second approximation for the equation of motion is 

F = /-^ (r, r) + h^2} ''^)- (9) 

The same procedure can be repeated iteratively for all higher orders. If the sequence of 
terms 6|„}(r, r) converges to a b{r,f) — in some appropriate sense, e.g., in a pointwise 
sense — then we arrive at a second order equation of motion, 

^ f(r,f) + P(r,r). (10) 

Naturally, it is a tough problem is whether the sequence in question converges or not. 



2.2 Iteration of the solution 

The above iteration method suggests another one, which would not directly result in an 
equation of motion but in the motion r corresponding to initial values Xq of spacetime 
position and Uq of four- velocity only. 

Let the solution of the zcroth approximation (3) corresponding to initial values xq 
and Mo be denoted by 2/{o}- Taking its first and third derivatives, we establish the first 
approximation 

F = fir, r) + V (^9\ - ^-h^^ y^y (11) 

Let Uj^ij be its solution for initial values Xq and Uq. Taking its third derivative, we establish 
the second approximation, and so on; we obtain 

= P{r, f) + 77 [g^ - ^i^fe^^ yl^^^, (12) 



If the sequence of solutions y^n} converges to a y then we succeeded in finding a motion 
satisfying equahty (2) without the need for the initial value of acceleration. This motion 
is, therefore, a good candidate for the sought physical solution. 

We can observe that, although this second method provides solutions rather than 
the equation of motion, a corresponding self-force function b{r, f) can be read off from 
the solutions. Namely, the value of the corresponding h at any spacetime point xq and 
four- velocity value uq can be calculated as 

V{xo,uo)=f,y'{Q), (13) 

i.e., taking the third derivative of the solution y belonging to initial values xq and uq, at 
the initial proper time value. 

Naturally, here, too, convergence is a hard problem. 



2.3 Equation for the self-force function 



Instead of utilizing some physically encouraging iteration scheme, we also have the possi- 
bility to directly formulate the desired property of the self- force function. Indeed, comput- 
ing the third derivative from the expected equation of motion (10) and then substituting 
f by the rhs of (10), the result demands from the self-force 



(14) 



This is a partial differential equation (system) for the two- variable function (s) lP{r^f). 
Naturally, one has to keep in mind that, here, r and r are to be understood as independent 
variables. We can display this fact in a more self-explaining way, writing (14) as 



ljP{r,u)^fi\g>y,- 



U^Uk 




du'- 



(15) 



We have arrived at a first order partial differential equation. Its solution is expected 
to contain an arbitrary free function, which arbitrariness is to be fixed to obtain the 
sought self-force function uniquely. 

On physical grounds, one can impose some requirements, via which one can rule out 
the undesirable solutions. We are going to use two such conditions. 

To formulate the first one, let us draw attention to that the self-force depends on the 
external force, b{r,r) = bf{r,r). Then we wish to demand that, if the exterior action 
tends to zero, then the self-force must tend to zero, too. Specifically, we will consider the 
self- force function b^f for every < n < I, and will impose 



limfe.^(r,r) 







(16) 



in the pointwise sense. 

The other assumption we will utilize below is that, if the exterior field has a spacetime 
symmetry, then the self-force function has the same symmetry. Namely, if / is invariant 
under a Poincare transformation P (with the underlying Lorentz transformation L), i.e., 
L~^f{Pr, Lf) — f{r, r), then the same invariance must hold for bf, too. 



3 Applying the three approaches in the nonrelativis- 
tic regime: Constant field 



Let us consider a constant external electromagnetic field, which acts on the charged 
particle via the Lorentz force 

f = eE + ev X B = eE + eFv, (17) 

where E is the electric field three- vector, B is the magnetic axial vector field, F = (— Bx) 
is the corresponding antisymmetric three-tensor — now each assumed space and time 
independent — , and v is the velocity of the particle. With the shorthands 

E:=— E, B := — B, F := — F, (18) 
m m m 

we can simply write 

f = f (v) = E + Fv. (19) 

The electromagnetic nature of the field will not play any role here so the subsequent 
considerations will be applicable for any force of the form (19), including a constant 
gravitational attraction, for example. 

We will need some technical remarks regarding F. Its kernel is spanned by B, its 
range is the plane orthogonal to B, and, with P denoting the orthogonal projection onto 
this plane, we find 

F = PF = FP = PFP, (I - P)F = 0, F^ = -E^P, (20) 

where I stands for the three-identity tensor (and B is the magnitude of B). 

To keep the formulae shorter and more easily accessible, we first treat E = 0. In this 
case, the equation of motion without radiation is 

V = Fv. (21) 

The results for the general case E 7^ are presented in Sect. 3.4. 

3.1 Iteration of the radiation term 

It follows from the equation without radiation — the zeroth approximation (21) — that 
V = Fv = F^v. Accordingly, the first approximation of the radiation term is 

b{i} = L{i}V, L{i} = -TyS^P, (22) 

with which the first approximation of the equation of motion becomes v = ^F -|- L^ij^v. 
The radiation term derived from this equation equals 

b{2} = L{2}V, L{2} ^fj(F + L{i})'. (23) 
Repeating this again and again, we recognize the generic recursion formula 

b{n} = L{n}V, L{„} = ?7(f + L{„_l})^ (24) 



Since 



L{2} = fj{F + L{i})' = fj (-B^P - 2f)B^F + fB'p) , 



(25) 



every L{„} is a linear combination of F and P. 

Supposing that the sequence L{„j. converges to an L since b{„}(v) should converge 
for any v , we have 



L = 7y(F + L)2, 
where L is a linear combination of F and P, 

L = -aF - 7P. 

Consequently, 

-aF - 7P = 77 [-(1 - afB^ + 27(1 - a)F + 7^?] , 
which imposes the pair of equations 

a = 2777(1 — a), 7 = {fjBy^l — ay — fjj'^. 
The first equation excludes a — 1, and then another equivalent pair of equations is 



?77 = 



2(1 -a)' 



4(77^)2(l-a)^+(l-a)2-l = 0. 



(26) 



(27) 



(28) 



(29) 



(30) 



The latter condition is a quadratic equation for ^ := (1 — a)^ > 0, with the only non- 
negative root 



Q 



-1 + ^1 + 16(77^)2 



(31) 



Prom here, we derive a = 1 ± ^/g. 

We can observe that two solutions are found. If convergence holds then, naturally, 
only one of them is the limit. To select the correct one, we can consider the simple case 
when there is no external force. Then all the iteration terms are zero, and the iteration 
converges trivially. Hence, the coefficients must be zero for zero magnetic field: 



q; = and 7 = for S = 0. 



(32) 



This requirement is, actually, in accord with the physical condition already worded in 
Sect. 2.3 that, for vanishing external field, the radiation term should vanish. It is easy to 
find from (31) that ^ — >■ 1 as S — > 0, and then (32) resolves the ambiguity and provides 
the final result 



a 



7 



2?7 \VQ 



1 



(33) 



the latter obtained from (30). 



The resulting equation of motion, with the radiation term (27) containing the coeffi- 
cients (33), reads 



V = (1 — a)F — 7P V, i.e., v = Fv + b(v) with b(v) = —aF — 7P v. 

(34) 

It is informative to inspect the solution of this equation of motion, which is 

v(t) = (I - P)vo + e-^*e('-")^*Pvo (35) 

for initial velocity Vq at zero initial time. For nonzero B,0<1 — a<l and 7 > 
so radiation causes that the effect of the external magnetic field is decreased and the 
component of velocity perpendicular to the magnetic field tends to zero as time passes. 

A problem with the method of iterating the radiation term is that it is difficult to 
obtain conditions for the convergence of the iteration. Unfortunately, convergence does 
not hold necessarily. Indeed, for fjB = 1 — i.e., for such special magnetic fields — 
convergence does not occur, since then L{2} = — 2F, so L{3j. = L^i}. Consequently, for 
all higher ns. 



L{2} {n is even), 
Ljij. (n is odd). 



(36) 



It is interesting, however, that the final result (33) with (31) does not exclude the case 
r]B — 1, and provides solution for rjB > 1 as well (which is expected to be outside the 
domain of convergence). 



3.2 Iteration of the solution 

Considering an initial value Vq at zero initial time, the zeroth equation (21) has the 
solution 

W{o}(t)=e^Vo. (37) 
Then W{o}(t) = — S^e^*Vo, so the first approximation satisfies the equation 

V = Fv - fjB^e^^o, (38) 
whose solution with initial value Vq is 

W{i}(t) = {l-fjBH)e^'Yo. (39) 

Then w{i}(t) = -[b^ (l — fjBH^ +2j7i?^Fje^*vo. The second approximation satisfies the 
equation 



v = Fv-fj[B^(l- fjBH) + 2^^^F 



Ff 



with the solution 



W{2}(t) 



l-fjB^(t- -fjBH ) - 2fj^BHF 



e^Vo. 



(40) 



(41) 



(42) 



At the nth step of this iteration scheme, we find 

W{n}(t) = [P{n}{t) + q{n}{t)F 

where and q^n} ^ire polynomials of t satisfying the recursive formulae 

P{n} = fj {p{n-\} - '2.B'^q{n-\} " B'^Pn-l) , ?{n} = V (?{n-l} + 2p{„_i} - B'^q{n-1}) (43) 

and the initial conditions p{„}(0) = 1, 9{n}(0) = 0. 

Let us suppose that p := \imnP{n} ^iid q := lim„g{„} exist and, moreover, that the 
limit procedure and difi^erentiation can be interchanged. Then we have for w := lim„ W{„} 

(44) 



w{t) = \p{t) + q{t)F 



— fj (j) — 2B^q — B^p^ , q — fj (q + 2p — B^q^ 



P 



together with the conditions p(0) = 1, g(0) — 0. With the notation s : = 
can be comprised as 



(45) 

p + iBq, (45) 



s - - 2i5J s-B^s^ 0. 
The roots of the corresponding characteristic polynomial are 



(46) 



— -iB±. 
2fj \ 



4:ff 



B 

1 — . 

V 



(47) 



The emerging ambiguity can be resolved as previously: if the magnetic force is sent to 
zero then p|„}(t) = 1 for all n, therefore, we choose the root that vanishes for ^ so 
as to obtain the solution p{t) = 1, i.e., the free motion (no self-force). 

Therefore, if (—7) and a are the real and the imaginary part of the root, respectively, 
then the solution of our system of differential equations is 



/X -^f / V /V ^.sm(aBt) 

p{t) = e-^* cos{aBt), q{t) = -e~^* — — '-, 



(48) 



where 



a — 2777(1 — a), 7 + 777^ = fjB^{l — aY- 



(49) 



Comparing the results with the outcomes of the approach of the iteration of the 
radiation term reveals that equations (49) coincide with (29). Furthermore, utilizing 
(20), one can demonstrate that (44) is the same as (35). At last, it is straightforward 
to find that the self-force function corresponding to (44) — see (13) — proves to be the 
same as (34). On the other side, one can also observe that, though leading to the same 
result, the two iterations themselves are different: the solutions of the iterated equations 
V = -|- L{„}^v do not equal the functions (42). 

Naturally, convergence is a nontrivial question in this approach, too. Nevertheless, it 
is interesting that, with this method, there is no evidence for excluding the case 77^ = 1. 



3.3 Equation for the self- force function 



We can start with ruling out the space and time dependence of 13, based on the require- 
ment of 2.3 that a spacetime symmetry of the external force should be respected by the 
self-force, too. In the present case, the symmetry in question is spacetime translation 
invariance. Hence, the sought equation of motion is of the form 



v = Fv-Fb(v). 



(50) 



In addition, the exterior field is invariant for space inversion, — F(— v) = F(v), and thus 
— b(— v) = b(v) is required, too. 

According to our third method, described in Sect. 2.3, the differential equation b must 
obey 



b(v) =f]\F + b'(v)l [Fv + b(v) 



(51) 



with ' denoting the derivative map. 

Let us assume that we can expand b in a series. Because of the space inversion 
symmetry, the even powers are zero, so 



b(v) = Liv + L3(v, V, v) + L5(v, v, v, v, v) + 



(52) 



where Li is a linear map, L3 is a symmetric trilinear map etc.; keep in mind that they 
depend on F. Using the notation Ls^^v^j := L3(v, v, v) etc., we obtain 

Liv + L3(v3) +L5(v^ 



F + Li + 3Ls(v^ ■ ) + 5L5(v^ • ) + • • 
from which it follows, order by order, that 

Li = 7?(f + Li)', 



Fv + Li V + L3 (v^) + L5 (v^) + 



L3(v^) = fj[3Ls{v^ (f + Li)v) + (f + Li)L3(v3) 



= 77[3L3(v^ L3(v2)) + 5L5(v^ (f + Li)v) + (f + Li)L5(v^) 



(53) 



(54) 



(55) 



(56) 



etc. Multiplying (54) from the left and from the right by F-I-Li, we find that LiF = FLi. 
Then it is a simple algebraic fact that 



Li = -q;F-7P + /3(I-P), 



(57) 



where a, 7 and P are scalar coefficients depending on F. Applying the first condition put 
in Sect. 2.3, 7 and j3 must tend to zero if F tends to zero. 
Having (57), (54) yields 

-aF - 7P + ^(I - P) = 77 [-(1 - a)'^B^P + 7^P + /3^{I - P) - 27(1 - a)F] , (58) 
which tells 

a = 2777(1 - a), -f^fi{l-afB^-fi-f^ (59) 



as well as /3 = fj^'^, according to which (3 is either zero or equals 4. The second possibility 
is excluded by the fact that (3 must be zero for zero external field so 

^ = 0. (60) 

Now we can recognize that (59) equals (29). As a result, we have (33) for the explicit 
expression of a and 7. 

As it is proved in the Appendix, (55) yields L3 = 0, and, similarly, all the higher 
order terms are found to be zero. Hence, the self-force is 

b(v) = -(q;F + 7P)v, (61) 

and the equation of motion is (34). 

Note that, in this simple case, the two conditions given at the end of Sect. 2.3 suffice to 
determine the self-force function completely, and we obtain the result derived by iterations 
without any problem on convergence. 



3.4 Nonzero electric and magnetic field 

As anticipated in Sect. 3, now we turn towards the case of a nonzero constant electric 
field in addition to the nonzero constant magnetic field. 

It is beneficial to decompose the electric field into components parallel to and orthog- 
onal to the magnetic field, respectively, and to observe that this decomposition can be 
written in the form 

E= (I-P)E + F (^--i^FE^ . (62) 

This induces a decomposition of the equation of motion without self-force [i.e., containing 
only the external Lorentz force (19)]. The B-parallel component, 

[(I-P)v]- = (I-P)E, (63) 

governs only the B-parallel component of v, while the B-orthogonal part can be written 
as 

li = Fu with u Pv - iFE, u ± B, (64) 

and determines the time evolution of the B-orthogonal component of v. 

Now we add the self-force term. Both iteration methods, namely, that of the radiation 
term and that of the solution, can be found to preserve this decomposition, where the 
B-orthogonal part can actually be treated the same way for u as we proceeded in the 
E = case for v. Both approaches provide 

b(v) = -aPE + T^FE - \aF + 7PI v (65) 

L J 

for the self-force and 

v(t) = (I - P)Et + (I - P)vo + iFE + e-^*e(i-")^* (pvq - -^Fe) (66) 

for the solution, after putting the two decomposed parts together. The found self-force 
function proves to be a solution of the partial differential equation of the third approach 
as well, and satisfies the two additional requirements — vanishing for vanishing external 
field and symmetry preservation. 



4 Applying the three approaches in the nonrelativis- 
tic regime: Elastic force 



As our second physical system considered as example for the three proposed methods, 
we next investigate the one-dimensional nonrelativistic motion due to a harmonic elastic 
force. Without radiation back-reaction, the equation is 

X — —uP'x^ (67) 

where a; is a non- negative constant. 



4.1 Iteration of the radiation term 

It follows from the equation without radiation (the zeroth approximation) that the first 
approximation of the radiation term f\x equals —fiuP'x. Accordingly, the first approxi- 
mation of the radiation self-force term is 

6{i}(a;, i) = —fjijj'^x, (68) 

and the first approximate equation of motion becomes x — —uj^x + b^iy{x, x). Computing 
fj'x from this equation and then replacing x with —uj^x + b{iy{x, x), we obtain the second 
approximation 

6|2}(x, x) = ffu'^x — fjoj'^ (l — ffoj'^^ X (69) 

and the second equation x = —oJ^x + 6{2}(x, i). 

It is straightforward then that the nth approximation is of the form 

(x, x) = anUj'^x - -fnX, (70) 

where q;„ and 777^ are functions (polynomials) of 77^0;^ and satisfy the recursive formulae 

q;„+i = 777„(1 - «„), 777n+i = 77^a;^(l - «„) - (777n)^. (71) 

Supposing that the sequence converges to an b, i.e., the limits a :— lim„Q;„ and 
7 := lim„7„ exist, then 

b{x, x) — auj^x — ^x (72) 

and 

« = T^, 777(1 + ^7)' = 77 V. (73) 
1 + 777 

We can see that 7 > and < a < 1 are necessary for a; 7^ 0. 

The roots of the second equation in (73) give 777 as a function of 77^0;^. Transforming 
that equation into the form (777 + 2/3)^ - (1/3) (777 + 2/3) - (2/27 + 77^0;^) = 0, we can 
find its solutions (roots) by the Cardano formula. Naturally, we are interested in the real 
roots. We find that only one root is real; with the notation 



the real root is 

2 

f)7 = + + (75) 

Hence, the hmiting equation of motion, with the radiation term (72), reads 

X — —(1 — a)uP'x — ^x. (76) 

According to this equation, radiaton causes that the effect of the elastic force is decreased 
and the motion is damped. 

It is a problem that it is difficult to obtain conditions for the convergence of the 
iteration. Unfortunately, convergence does not hold necessarily. Indeed, if 77a; = 1 then 
the sequence does not converge because then 6(2} (a;, x) = , and thus 6{3} = 
Consequently, 

t /^{2} (n is even), 
I 0{i} [n IS odd). 

It is interesting, however, that, with (75), (73) provides a solution for all values of f\LO. 
4.2 Iteration of the solution 

Considering some initial values Xq and for position and velocity, the zeroth equation 
(67) has the solution 

y{o}(i)=-2oe'"* + ^oe~'"* = -2oe^"* + c.c. with := ^ (2^0 - ^) • (78) 
Since j/{o}(^) — —izou^e"^^ + c.c, the equation for the first iterated solution reads 

X = -oj'^x + fj (-i^o^^V'^* + c. c.) , (79) 

whose solution with the chosen initial values is y{i}{t) — [— {1/2) fjouH + Zo]e"^^ + c. c. . 
Then we find that the nth solution is of the form 

?/{n}(i) =Pn(i)e'"* + c.c. (80) 

where p„(t) is a polynomial of t of nth degree; then 

y{n} = (Pn + 'iiupn - V)e"^* + C C. (81) 

and 

y'in} = (Pn + ^iojpn - ?>u'^pn - ioJ^PnW''^ + C. C. , (82) 



and the resulting recursive formula (arising both from the coefficient of e"^* and from that 
of e-'^*) is 

Thus, supposing convergence (and that differentiation can be interchanged with taking 
the limit), we find for the limit p := lim„p„ 

O O /*** O Q\ /\ 

p + 2iujp — u! p = —ui p + fj{P + Siuip — 3ui p — iui p). (84) 
This is a hnear differential equation of third order. Its solutions are of the form e"^* with 

{X + iuf + u'^ = f}{X + iuf. (85) 
Putting \ + iu) = 11 + iu, where p and u are real, we can rewrite this as 

+ 2ipLu - + a;^ = f}{n^ + 3i//V - Spiu'^ - iV^), (86) 
which is equivalent to the pair of real equations 

i/2^3//2-^, 8{fifj,f -8{wf + 2{fjfx)- {fjcof ^0. (87) 
V 

With 7 := 2/x, the latter reduces to the second equation of (73). In parallel, putting 
v'^ —: (1 — a)uj^ — 7^/4 (i.e., defining a in this way), we find fij{fij + 1) = (1 — a)u^, 
which, together with the second equation of (73), results in the first one. 

As a consequence, whenever convergence holds, the iteration of solutions gives the 
same solutions and same self-force function as the iteration of the radiation term. 

4.3 Equation for the self-force function 

According to the third method (Sect. 2.3), the radiation self-force function b{x,x) in the 
anticipated equation of motion 

X ^ -uj'^x + b(x,x) (88) 
is to satisfy the quasi-linear partial differential equation 



b{x, v) — rj 



-uh + ^^^v + ^^(-u^x + b) 
ox OV ^ ' 



(89) 



Note that, although not denoted explicitly, the sought h also depends on a;, i.e., on the 
external force. 

The characteristic ordinary differential equation corresponding to (89) reads 

dx di) 9 - d6 1 - 9 , , 

— = —^-u?x + h, —^-h + u?v. 90 

d^ 77 

This is a simple linear differential equation whose characteristic roots A fulfill the equation 
77A^ — A^ — a;^ = 0. For 777 := 77A — 1, this is exactly the second equation of (73). Thus, 
with the notation (74), we have the roots 



^As = - + p+ + p- 



(92) 



With the three roots, the solutions of equation (90) are of the form 

x{0 = E Qe^^^ v{0 = E A,c,e^'^ 6(0 = E A?c,e^'« + co'xiO- (93) 

i=l i=l i=l 

According to the method of characteristics, via ehminating the auxihary variable ^, one 
obtains 6 as a function of x and v. Now we use the condition of Sect. 2.3 that b must be 
zero for zero external force, which in the present case means a; = 0. Since Ai and A2 are 
zero for a; = and A3 is not zero, we deduce from the last equality above that C3 = is 
necessary. 

Then ^ can be eliminated easily; with the notations Zi :— Cie^^^ and Z2 '■— C2e'^^^, we 
have 

X ^ zi + Z2, V ^ XiZi + X2Z2, b = X^zi + \\z2 + uP'x. (94) 

Here, the first two equations enable one to express zi and Z2 as a linear function of x and 
V. Substituting them into the third equation provides 6 as a linear function of x and v: 
b{x, v) = {uj^ - XiX2)x + (Ai + X2)v. 

A convenient way to proceed is to write the coefficients in another form: 

b{x,v) = au^x — jv. (95) 

Substituting this form into (89), we obtain (73) for a and 7. Consequently, we can derive 
the same conclusions as in Sect. 4.1. 

We can see that in this case, too, the condition that the self-force must be zero if 
the external force is zero suffices to determine the self-force function completely, and 
we recover the result derived by the previous two methods without any problem on 
convergence. 

5 Discussion 

We have suggested three methods to determine a second order equation of motion whose 
solutions satisfy the Lorentz-Dirac equality and, at the same time, are physically ac- 
ceptable. Our aim was to give back-reaction as a function of spacetime position and 
velocity. 

In the nonrelativistic approximation of the Lorentz-Dirac equality, we could exactly 
calculate this force function for two systems: a constant external electromagnetic field 
and a one-dimensional elastic external force. The three methods turned out to lead to 
the same back-reaction force function. According to the resulting physical picture, for 
both systems, radiation back-reaction has two manifestations: inducing a damping linear 
in velocity, and decreasing the strength of the external force. 

The latter effect could also allow the — quantum field theory motivated — interpre- 
tation that back-reaction causes a positive renormalization of the mass of the particle. 
However, for the constant external electromagnetic field, this renormalization turns out 
to differ from a simple scalar multiplying of the mass. Rather, renormalization must be 
a direction dependent, tensorial multiplication. 



This system also proves to show a hmitation of the criterion by Dirac and Haag [6, 7], 
which would choose that solution for initial position and velocity for which acceleration 
tends to zero for asymptotically large times. In fact, this system is found to decouple into 
two independent subsystems, one parallel to and the other orthogonal to the magnetic 
field. In the former subsystem, acceleration remains time independent and corresponds 
to unrenormalized mass, while damping and renormalized mass (or renormalized external 
force) emerges in the latter subsystem. 

It is therefore an important open question on what grounds the decrease of the ex- 
ternal force can be interpreted as mass renormalization, and whether in other systems 
this renormalization is not only direction dependent but, for example, also (spacetime) 
position and velocity dependent (as suggested by some preliminary considerations not 
detailed here). 

The three methods we proposed and investigated led to the same result for the two 
systems considered. Some differences among the three approaches were found, though. 
First, in the two iteration approaches, not all coefficients of the model ensured the exis- 
tence of a fixed point. Second, convergence to the found fixed points remained a tough 
open mathematical problem. On the other side, in the third approach, the problematic 
aspect is that there is an encoded amount of ambiguity in the partial differential equation 
to solve. This ambiguity was easy to rule out for the two systems we considered but may 
be a harder task for other systems. 

Further study is needed, accordingly, about each method separately and also about 
some possible connections among them. 



Appendix 

We start with deriving a simple algebraic fact. Let T be a symmetric trilinear map. Then 
(with the simplifying notation used earlier) 



t((w + v)^) = t(w^) + 3T(w^v) + 3T(w,v2) + t(v^). 



(96) 



Thus, if t(( • )^) = then 3T(w2, v) + 3T(w, v^) = for all v and w. Here, for a fixed 
V, the first term is bilinear, the second term is linear in w; their sum can be zero only if 
both are zero. As a consequence, if T^v^j = for all v then T(w^, v) = for all w and 
V. Further, for a fixed v, we have 



T((wl + W2)^v) = t(w?,v) +2T(wi,W2,v) + t(w^,v), 
which shows that, if T^w^, = for all w and v, then 

T(wi, W2, v) = for all Wi, W2, v. 

Hence, we have the final result: if T(v^) = for all v then (98) holds. 
The previous consideration and (55) yield that 

I-57(F + Li) L3(viri, W2, v) = ST^La^wi, W2, (F + Li)v) 

for all Wi, W2 and v. 

For fixed wi and W2, W := L3(wi, W2, ■ ) is a linear map satisfying 



(97) 



(98) 



(99) 



I-f}{F + Li) W = 3fjW{F + Li), 



(100) 



which imphes 



W 



I - 3r^(F + Li) = (F + Li)?7W, (101) 



too. Recall that F + Li = (1 — a)F — 7P. Then it is a simple fact that the hnear maps 
multiplying W on the left hand side in (100) and (101), respectively, arc nondcgcncrate. 
Therefore, applying I — P from the right to (100) and from the left to (101), we find 

WP = W = PW. (102) 

As a consequence, W commutes with F, too, so it must be of the form 

W = ^F + CP, (103) 

and either (100) or (101) implies 

(1 + 4777)W = 477(1 - a)WF. (104) 

This equation and the previous one result in 

(1 + 4777)/? = 477(1 -q;)C, (1 + 4777)C = -477(1 -q;)^^^ (105) 

implying 5^/3^ = — C^, which is possible only if /3 = and C = 0) i-G-) W = 0. 
Since we had W = L3(wi, W2, • ) for arbitrary wi and W2, we arrive at 

L3 = 0. (106) 
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